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Abstract 

We investigate D- dimensional gravitational model with curvature-quadratic and curvature-quartic cor- 
rection terms: R + R 2 + R 4 . It is assumed that the corresponding higher dimensional spacetime manifold 
undergos a spontaneous compactification to a manifold with warped product structure. Special attention 
is paid to the stability of the extra-dimensional factor space for a model with critical dimension D — 8. It 
is shown that for certain parameter regions the model allows for a freezing stabilization of this space. The 
effective four-dimensional cosmological constant is negative and the external four-dimensional spacetime is 
w**) ' asymptotically AdS. 

o' 

1 Introduction 

Mh. Multidimensionality of our Universe is one of the most intriguing assumption in modern physics. It follows from 

theories which unify different fundamental interactions with gravity, such as M or string theory £Q, and which 
have their most consistent formulation in spacetimes with more than four dimensions. Thus, multidimensional 
cosmological models have received a great deal of attention over the last years. 

Stabilization of additional dimensions near their present day values (dilaton/geometrical moduli stabiliza- 
tion) is one of the main problems for any multidimensional theory because a dynamical behavior of the internal 
spaces results in a variation of the fundamental physical constants. Observations show that internal spaces 
should be static or nearly static at least from the time of recombination (in some papers arguments are given 
in favor of the assumption that variations of the fundamental constants are absent from the time of primordial 
nucleosynthesis [2]). In other words, from this time the compactification scale of the internal space should either 
be stabilized and trapped at the minimum of some effective potential, or it should be slowly varying (similar to 
the slowly varying cosmological constant in the quintessence scenario). In both cases, small fluctuations over 
stabilized or slowly varying compactification scales (conformal scales/geometrical moduli) are possible. 

Stabilization of extra dimensions (moduli stabilization) in models with large extra dimensions (ADD-type 
models) has been considered in a number of papers (see e.g., Refs. [H]-^] )• In the corresponding approaches, a 
product topology of the (4 + D')— dimensional bulk spacetime was constructed from Einstein spaces with scale 
(warp) factors depending only on the coordinates of the external 4— dimensional component. As a consequence, 
the conformal excitations had the form of massive scalar fields living in the external spacetime. Within the 
framework of multidimensional cosmological models (MCM) such excitations were investigated in m]-E3 where 
they were called gravitational excitons. Later, since the ADD compactification approach these geometrical 
moduli excitations are known as radions [<2|S]. 

Most of the aforementioned papers are devoted to the stabilization of large extra dimensions in theories with 
a linear multidimensional gravitational action. String theory suggests that the usual linear Einstein-Hilbert 
action should be extended with higher order nonlinear curvature terms. In our papers ^7]-[2D] we considered a 
simplified model with multidimensional Lagrangian of the form L = f(R), where f(R) is an arbitrary smooth 
function of the scalar curvature. Without connection to stabilization of the extra-dimensions, such models 
(4— dimensional as well as multidimensional ones) were considered e.g. in Refs. |14j-jlri|. There, it was shown 
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that the nonlinear models are equivalent to models with linear gravitational action plus a minimally coupled 
scalar field with self-interaction potential. Similar approach was elaborated in Refs. |2J where the main 
attention was paid to a possibility of the late time acceleration of the Universe due to the nonlinearity of the 
model. In our papers |17 | -|2(J | . we advanced the equivalence between the nonlinear models and the linear ones 
with a minimally coupled scalar field towards investigating the stabilization problem for extra dimensions 1 ). 
Particular attention was paid to models where usual linear curvature term R was supplemented with either 
R 2 or R 4 or R^ 1 nonlinear terms. All of these models can be investigated analytically. It was shown that for 
certain parameter ranges, the extra dimensions are stabilized. In the present paper we extend this consideration 
to a model with Lagrangian of the type f(R) = R + R 2 + R A - Such simple generalization enriches considerably 
the qualitative behavior of the model because it results in either one-branch or three-branch models. For each 
of these branches the stability analysis should be performed separately. Our paper is the first one among two 
papers devoted to this model. Here, we investigate the one-branch model postponing the three-branch model for 
our forthcoming paper. This model can be solved analytically. However, this analysis is very cumbersome for 
an arbitrary number of dimensions D. So, we restrict our consideration to a critical dimension D = 8. Critical 
dimension is defined by doubled degree of the scalar curvature polynomial f(R), i.e. in our case D = 2 x 4 
(see |2S] and Appendix A below). We show that for certain parameter regions the model allows for a freezing 
stabilization of the internal space. Here, the stabilization takes place in negative minimum of the effective 
potential. Thus the effective four-dimensional cosmological constant is also negative and the homogeneous and 
isotropic external four-dimensional spacetime is asymptotically AdS. 

The paper is structured as follows. In section we present a brief technical outline of the transformation 
from a non-Einsteinian purely gravitational model with general scalar curvature nonlinearity of the type f(R) 
to an equivalent curvature-linear model with additional nonlinearity carrying scalar field. Afterwards, we derive 
criteria which ensure the existence of at least one minimum for the effective potential of the internal space scale 
factors (volume moduli). These criteria are then used in Sec. 01 to obtain the regions in parameter space which 
allow for a freezing stabilization of the scale factors in the model with scalar curvature nonlincaritics of the 
type f(R) — R + R 2 + R A . A brief discussion of the obtained results is presented in the concluding Sec. 0] 
In appendix lAl it is shown that D = 2 x deg^(/) is the critical dimension of the considered models. Finally, 
graphical visualization of the effective potential with a global minimum is given in appendix [5] 

2 General setup 

We consider a D = (4 + D')— dimensional nonlinear pure gravitational theory with action functional 
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where f(R\ is an arbitrary smooth function with mass dimension 0(m 2 ) (m has the unit of mass) of a scalar 
curvature R = R[g] constructed from the D~ dimensional metric g a b {a,b = 1, . . . , D). D' is the number of extra 
dimensions and k 2 d denotes the D— dimensional gravitational constant which is connected with the fundamental 
mass scale M*u+d') and the surface area Sd-i = 2tt^ d ^ 1 ^/ 2 /T[(D — l)/2] of a unit sphere in D — 1 dimensions 
by the relation (2HI HZ] 

K D = 2S D -i/M^ i+D ,y (2.2) 

Before we endow the metric of the pure gravity theory l|2.1|l with explicit structure, we recall that this 
R— nonlinear theory is equivalent to a theory which is linear in another scalar curvature R but which con- 
tains an additional self- interacting scalar field. According to standard techniques |141 ITS] , the corresponding 
R— linear theory has the action functional: 

s = ^T [ d D xV\g~\ [R[g] ~ g ah <$>^. h - 217(0] , (2.3) 

where 



S =^T I d D x^/\f\f(R), (2.1) 



% K D JM 



m- &-■»>». -W§^' (2 - 4) 

and where the self-interaction potential U((j>) of the scalar field <f> is given by 

U{4>) = \{ir D/{D ~ 2) '[Rf '- f] , (2.5) 

= le-^ \ R(ct>)e A * - f (R(<f>))] , B := — = (2.6) 
2 L ^ V yV ' n v /(£>-2)(D-l) 



) A different approach to the problem of the extra dimension stabilization was proposed in Ref. 1221 . This method can be 
applied to Lagrangians containing high-order curvature invariants. In the case of the models with Lagrangians L = f(R) both of 
these approaches result in the same conclusions. 



This scalar field cj> carries the nonlinearity degrees of freedom in R of the original theory, and for brevity we 
call it the nonlinearity field. The metrics g ab , Sab of the two theories (|2.1[) and (|2.3|) are conformally connected 
by the relations 

ga b = n 2 -g ab =[f(R)f {D - 2) g ab . (2.7) 

As next, we assume that the D-dimensional bulk space-time M undergoes a spontaneous compactification 
to a warped product manifold 

M = M x M x x . . . x M n (2.8) 

with metric 

n 

g = g ab (X)dX a <g> dX b = g^ + ^ e 2/5l(x) 5 (l) . (2.9) 



; = 1 



The coordinates on the (D a = rf + 1)— dimensional manifold M (usually interpreted as our observable (D = 
4)— dimensional Universe) are denoted by x and the corresponding metric by 

g(°) = g$(x)dx^ ® dx v . (2.10) 

For simplicity, we choose the internal factor manifolds Mi as di— dimensional Einstein spaces with metrics 
9^ — 9rn i n i {Vi)dy r [ li (g) dy™ 1 , so that the relations 



Rr, 

and 



9® 



= \ i g$. n ., m i) n i = l,...,di (2.11) 



R 



9® 



= Xdi = Rt (2.12) 



hold. The specific metric ansatz (J2.9II leads to a scalar curvature R which depends only on the coordinates x of 
the external space: R[g] = R(x). Correspondingly, also the nonlinearity field <p depends on x only: <j) = (j)(x). 

Passing from the R— nonlinear theory l|2.1() to the equivalent R— linear theory l|2.3|l the metric l|2.9|l undergoes 
the conformal transformation g i— > g [see relation (|2.7|l ] 

n 

g = n 2 g= (e A *) 2/iD - 2) g := g™ + ]T e^^g^ (2.13) 

with 

5 (o) ;= ( ^/M - g (0) t ^ := pi + _A_ . (2 . 14) 

The main subject of our subsequent considerations will be the stabilization of the internal space components. 
A strong argument in favor of stabilized or almost stabilized internal space scale factors f3 % (x) , at the present 
evolution stage of the Universe, is given by the intimate relation between variations of these scale factors and 
those of the fine-structure constant a |28j . The strong restrictions on a— variations in the currently observable 
part of the Universe [2H imply a correspondingly strong restriction on these scale factor variations [2U- For 
this reason, we will concentrate below on the derivation of criteria which will ensure a freezing stabilization of 
the scale factors. Extending earlier discussions of models with R 2 , i? 4 and R~ x scalar curvature nonlinearities 
^7]-|2ni w e will investigate here models of the nonlinearity type R 2 + R 4 . 

In Ref. J3I it was shown that for models with a warped product structure l|2.9|l of the bulk spacetime M 
and a minimally coupled scalar field living on this spacetime, the stabilization of the internal space components 
requires a simultaneous freezing of the scalar field. Here we expect a similar situation with simultaneous freezing 
stabilization of the scale factors (3 l (x) and the nonlinearity field <fi(x). According to (|2.14|) . this will also imply 
a stabilization of the scale factors /3 4 (x) of the original nonlinear model. 

In general, the model will allow for several stable scale factor configurations (minima in the landscape over 
the space of volume moduli). We choose one of them 2 ), denote the corresponding scale factors as /3q, and work 
further on with the deviations 

P{x)=P{x)-Pl (2.15) 

as the dynamical fields. After dimensional reduction of the action functional l|2.3() we pass from the intermediate 
Brans-Dicke frame to the Einstein frame via a conformal transformation 

/ n \ -V(A>-2) 

^=^^=(11^1 C (2.16) 

2 ) Although the toy model ansatz 12.11 is highly oversimplified and far from a realistic model, we can roughly think of the chosen 
minimum, e.g., as that one which we expect to correspond to current evolution stage of our observable Universe. 



with respect to the scale factor deviations $ l (x) J7j-|2nij|23- As result we arrive at the following action 
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d D °, 



lff (0) l{ 



R 



S<°> 



Gifg^d^ dji - g^d^d^ - 2J7 e //} 



(2.17) 
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which contains the scale factor offsets /3q through the total internal space volume 

n n 

V D , = V T x vo = TT / d di yJ\g^\ x J[ e d ^ 

in the definition of the effective gravitational constant n 2 D of the dimcnsionally reduced theory 

4tt .. . ,2+u' 



K (D =4) = «I>/Vd' = &7T/MI 



Ml = 



—V D> M <4+D , 



(2.18) 



(2.19) 



Obviously, at the present evolution stage of the Universe, the internal space components should have a total 
volume which would yield a four-dimensional mass scale of order of the Planck mass M( 4 ) = Mpi . The tensor 

components of the midisuperspace metric (target space metric on Kj.) reads: Gy = diSij +didj/(Do — 2), where 
i,j = (1, . . . ,n), see J3U1I2]- The effective potential has the explicit form 



U eff 0,ct>)= [He 



di,3* 



1 



IT,* 



-2/3' 



uw 



(2.20) 



where we abbreviated 

&:= Ri&cp (-20q). (2.21) 

A freezing stabilization of the internal spaces will be achieved if the effective potential has at least one 
minimum with respect to the fields l [x). Assuming, without loss of generality, that one of the minima is 
located at j3 l = /3q =>■ (3 l — 0, we get the extremum condition: 



dU, 



eff 



d/3 l 



= 



Rt = 



/3=0 



Dn 






(2.22) 



From its structure (a constant on the l.h.s. and a dynamical function of 4>{x) on the r.h.s) it follows that a 
stabilization of the internal space scale factors can only occur when the nonlinearity field 4>{x) is stabilized as 
well. In our freezing scenario this will require a minimum with respect to 4>: 



dU(<j>) 







dU, 



eff 



0. 



(2.23) 



We arrived at a stabilization problem, some of whose general aspects have been analyzed already in Refs. jllj - 
|l,'i| and |17 | -|2U ) . For brevity we only summarize the corresponding essentials as they will be needed for more 
detailed discussions in the next sections. 

1. Eq. (|2.22l) implies that the scalar curvatures Ri and with them the compactification scales e^° [see relation 
(|2.21|l ] of the internal space components are finely tuned 



Ri Rj 



hj = 1, • 



(2.24) 



2. The masses of the normal mode excitations of the internal space scale factors (gravitational excitons/radions) 
and of the nonlinearity field <j) near the minimum position are given as |13j : 



= m* = --^L^o)=-2§>0, 



d 2 U{cj>) 



>0. 



(2.25) 
(2.26) 



00 



3. The value of the effective potential at the minimum plays the role of an effective 4D cosmological constant 
of the external (our) spacetime Mq: 



A e// := U ef f 



/3'=0, 

= 00 



A) -2 
D-2 



Ufa) 



Do-2 R z 
2 di 



(2.27) 



4. Relation l|2.27|l implies 

sign A e// = sign U(<f> ) = signi?, . (2.28) 

Together with condition l|2.25|l this shows that in a pure geometrical model stable configurations can only 
exist for internal spaces with negative curvature 3 ): 



Ri<0 (i = l,...,n). 



(2.29) 



Additionally, the effective cosmological constant ^- e ff as well as the minimum of the potential U(4>) should 
be negative too: 

A e //<0, C/(0o)<O. (2.30) 

Plugging the potential U(<f>) from Eq. (|2.6(l into the minimum conditions l|2.23|l . I|2.26|l yields with the help of 
d^R = Af'/f" the conditions 



dU 

~dlp 

d 2 U 



A 



{fr n/ ( n-2) h 



<t>0 



2{D-2) 

1 Ae {A-B)4 r^g + ( A _ B) =j 



= 0, h:= Df - 2Rf, 
1 



2(^-i) (/) r di?/l 



M0o) = o, 
>o. 



4>o 



where the last inequality can be reshaped into the suitable form 

f"d R h\ <f)o =f"[(D-2)f-2Rf"] (f>B >Q 

Furthermore, we find from Eq. H2.31fl 

D-2 



so that (|2.30() leads to the additional restriction 



2D 



(/') D - 2 ii(0o) 



R(M < o 



(2.31) 

(2.32) 

(2.33) 

(2.34) 
(2.35) 



at the extremum. 

In the next section we will analyze the internal space stabilization conditions (|2.24|l - H2.30JI and l|2.31(l - 
(|2.35|) on their compatibility with particular scalar curvature nonlinearity f(R). According to our definition 
(|2.4|) . we shall consider the positive branch 

f'(R)>0. (2.36) 

Although the negative /' < branch can be considered as well (see e.g. Refs. ^JElEOl)) we postpone this 
case for our future investigations. 

3 The R 2 + i? 4 -model 

In this section we analyze a model with curvature-quadratic and curvature-quartic correction terms of the type 

f(R) = R : + aR 2 + 7 i? 4 - 2A D . (3.1) 

We start our investigation for an arbitrary number of dimensions D. First of all, we should define the relation 
between the scalar curvature R and the nonlinearity field <fi. According to eq. (|2.4Jl we have: 



/' = e M = 1 + 2aR + 4 7 i? 3 



(3.2) 



3 ) Negative constant curvature spaces Mi are compact if they have a quotient structure: Mi = H di /Ti, where H di and I\ are 
hyperbolic spaces and their discrete isometry group, respectively. 



This equation can be rewritten in the form 



R a + —R 



27 " 47 



= 0, 



where 



X = e A ' p - 1 , -00 < <p < +00 -£=4- -1 < X < +00 . 
Eq. Il.'i.."ill has three solutions Ri } 2,3, where one or three of them are real valued. Let 

q := — , r := — X. 
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The sign of the discriminant 

Q := r 2 + q 3 

defines the number of real solutions (see, e.g., Ref. |32p: 



(3.3) 

(3.4) 

(3.5) 
(3.6) 



Q>0 
Q = Q 
Q<0 



3i?! = 0, 3i? 2i3 + 

SsRi = 0,* = 1,2,3, Ri=R 2 
QR\ = 0,z' = l,2,3. 



(3.7) 



It is most convenient to consider Ri — Ri(X) as solution family depending on the two additional parameters 
(a, 7). Physical scalar curvatures correspond to real solutions Ri(X). For Q > the single real solution R\ is 
given as 



Ri= r + Q 1/2 
The three real solutions Ri^^X) for Q < read 
Ri = S1 + S2, 
Ri 



1/3 



)l/2 



1/3 



(3.8) 



si, 



Ra 



-(-1 + iV3)si + -(-1 - iV3)s 2 = e ,: T Sl + e - 
-(-1 - iy/S)si + -(-1 + iV3)s 2 = e-^si + e^s 2 , 



where we can fix the Riemann sheet of Q 1 ' 2 by setting in the definitions of si, 2 



si 2 := 



rizlQI 1 / 2 



1/3 



In this paper we investigate the case of positive Q(<j>) that is equivalent to the condition 

Q{4>) > => sign a = sign7 . 



(3.9) 



(3.10) 



(3.11) 



The case sign a =^ sign 7 that corresponds to different signatures of the discriminant Q will be considered in 
our forthcoming paper. 

To define the conditions for minima of the effective potential U e ff , first we obtain the extremum positions 
of the potential U(4>). The extremum condition l|2.31|l for our particular model (|3.1|l reads: 



^(0)i7 I y - 4) + R 2 m a ( - - 2) + R m ( £ - I I - £A D = , 



(3.12) 



where subscript 1 indicates that we seek the extremum positions for the solution (|3.8|) . Eq. (|3.12l) clearly shows 
that D = 8 is the critical dimension for the model (|3.1|l in full agreement with the result of the Appendix (see 
(IA.4|l ). In what follows, we investigate this critical case. For D — 8 eq. (|3.12|1 is reduced to a quadratic one 



■^(0)1 + g - -R(o)i — ; A 8 = A D=S 



with the following two roots: 



«=-iW(i) + ^ 



(3.13) 



(3.14) 



These roots are real if parameters a and As satisfy the following condition: 



3 \ 2 4 



, , +-A 8 >0. (3.15) 

4a / a 

If sign (a) = sign(Ag), then condition (|3.15|) is automatically executed, else 

|A 8 | < -£- , sign (a) + sign (A 8 ) . (3.16) 

o4|a| 

To insure that roots l|3.14|l correspond to a minimum value of U(4>), they should satisfy the condition i|2.33|) : 

/" [(D - 2)/' - 2i?/"] 0o > <=► /"[3 + 4ai?] 0o > , (3.17) 

where 

f" = 2a + 12~/R 2 . (3.18) 

Because for Q > eq. I|3.8|l is the single real solution of the cubic eq. 13.3fl . then R = Ri((f>) is a monotonic 
function of <\> Thus, the derivative d$R\ = Af'/f" does not change its sign. Keeping in mind that we consider 
the /' > branch, the function Ri(<f>) is a monotone increasing one for f" > 0. As apparent form eq. (|3.8|) . 
for increasing Ri we should take 7 > 0. In a similar manner, the function R\ (cf>) is a monotone decreasing one 
for /",7 < 0. Thus, for the minimum position .R(o)i) inequality (|3.17|) leads to the following conditions (we 
remind that according to eq. I|2.35|l the minimum position i£(o)i should be negative and according to eq. (|3.11J) 
sign a = sign 7): 



I. /" > 7,o>0 : 

II. /", 7 ,«<0 



3 + 4a J Rf 0)1 >0^|flf 0)1 |< — . (3.19) 

3 + 4aRf m < <=► -\Rf Q)1 \ > JL . (3.20) 

Obviously, inequality (|3.2UI) is impossible and we arrive to the conclusion that the minimum of the effective 
potential U e ff is absent if sign a = sign 7 = — 1. 

Additionally, it can be easily seen that in the case 

sign a = sign 7 = signA/) = +1 (3-21) 

the effective potential U e ff has no minima also. This statement follows from the form of the potential U(4>) for 
the model (|3.1() . According to eq. I|2.5|l . U{<f>) reads: 

U(4>) - (l/2)e- s ^ (aR 2 + 3 7 i? 4 + 2Ad) . (3.22) 

Thus, this potential is always positive for parameters satisfying (|3.21|) and we arrive to the contradiction with 
the minimum condition l|2.30|l . Therefore, the investigation carried above indicates that the internal space stable 
compactification is possible only if the parameters satisfy the following sign relation: 



a >0, 7 >0,A 8 <0. (3.23) 

Let us investigate this case in more detail. For this choice of signs of the parameters, it can be easily seen that 
Z/qxj from eq. I|3.14|l satisfy the condition (|2.35|) : R( Q \\ 



both extremum values -R[ H ^ TOln ec l- l|3.14|l satisfy the condition (|2.35|) : -R[ ~ii < ®- However, the expression 



?(±) j.>i«5(±)a_ 



1 , /9 



?(±) 



:i 



= 1 + 2aR^ ) [ + byRfig = -- ± \j 7 - 16a|A 8 | - 4 7 R™ (3.24) 

shows that only R( \\ can belong to /' > branch. To make /' I R( -* ) positive, parameter 7 should satisfy 
the condition 

~2 



7< 



! ;,i 16a|A 8 



-"•(o)i 



(3.25) 



4 
As apparent from this equation, parameter 7 remains positive if A 8 belongs to the interval 



\,G(-i-,0) . (3.L>.i) 



For this values of Ag, the condition ((3.16|) is automatically satisfied. We also note, that for positive a and 
negative iwgvl the condition (|3.19|) is also satisfied. Taking into account the interval 1(3.26(1 . the corresponding 
allowed interval for 7 reads 4 ) 

7 e (0, _ l {+) 3 1 • (3.27) 

^(0)1 

Thus, for any positive value of a, Eqs. ((3.26(1 and 1(3.27(1 define allowed intervals for parameters A 8 and 7 
which ensure the existence of a global minimum of the effective potential U e ff . Here, we arrive to the required 
stable compactification of the internal space. The position of the minimum (/?g, 4>q) and its value can be easily 
found (via the root R/q^) with the help of Eq. ((3.2|) and corresponding Eqs. from section 3. The Fig.l - Fig. 2 
(see Appendix B ) demonstrate such minimum for a particular choice of the parameters: a = 1, 7 = 1, Ag = — O.f . 
To conclude this section, we want to note that limit Ag — > corresponds to R( Q)1 — > which results in the 
decompactification of the internal space /3 1 , — » 00. 

4 Conclusions 

In our paper we analyze the model with curvature-quadratic and curvature-quartic correction terms of the type 
((3.1(1 and show that the stable compactification of the internal space takes place for the sign relation 13.23(1 . 
Moreover, the parameters of the model should belong to the allowed intervals (regions of stability) ((3.26(1 and 
((3.27(1 . The former one can be rewritten in the form 

A 8 = £, £e(-l,0). (4.1) 

Thus, for the root -R^u anc ^ parameter 7 we obtain respectively 

R$i = l> ^i(-3+v / 9T8e)<0 (4.2) 



and 



7=2rw. Ce(0,l). (4.3) 



Eq. (H2J shows that R^l G (-^,0). 

It is of interest to estimate the masses of the gravitational excitons 1(2.25(1 and of the nonlinearity field (j> 
((2.26(1 as well as the effective cosmological constant 1(2.27(1 . From Eqs. ((4.1(1 - 1(4.3(1 follows that -fwoli ~ ^ 8 ~ 
—or 1 , 7 ~ a 3 =^- f'(4>o) ~ 0(1) , f"((j>o) ~ a , U(<po) ~ —a~ l . Then, the corresponding estimates read: 

-A e// ~m\~m\~ a' 1 . (4.4) 

From other hand (see Eqs. 1(2.21(1 and ((2. 25( 1 1 

!7(to) ~ exp(-2$) = 6^ x . (4.5) 

So, if the scale factor of the stabilized internal space is of the order of the Fermi length: 6( )i ~ Lp ~ 10~ 17 cm, 
then a ~ L|, and for the effective cosmological constant and masses we obtain: — A e y/ ~ m^ ~ m^ ~ ITeV 2 . 

In the present paper the analysis of the internal space stable compactification was performed in the case 
Q(4>) > =>• signa = sign7. In our forthcoming paper we extend this investigation to the case of negative Q((jJ) 
where the function R((f>) has three real-valued branches. 
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A On critical dimensions in nonlinear models 

The existence of a critical dimension (in our case D = 8) i s a rather general feature of gravitational theories 
with polynomial scalar curvature terms (see, e.g., Refs. |221I^I2S1)- Following our paper |20], it can be easily 
demonstrated for a model with curvature nonlinearity of the type 

N 

f(R) = J2<*kR k (A.l) 

for which the ansatz 

N 

e A<t> = f = '£ka k R k - 1 (A.2) 

fe=0 

leads, similar like (|2.5() . to a potential 



1 
U{4>) = ~ {f')- D,{D - 2) J2(k - l)a k R k . (A.3) 



fc=0 



The condition of extremum (|2.31|) for this potential reads: 



JV 



Df - 2Rf = £ (-D - 2k) a k R k = . (A.4) 

fe=0 

Thus, at the critical dimension D = 2N the degree of this equation is reduced from N to TV — 1. In this case 
the search of extrema is considerably simplified. 

In the limit (f> — > +oo the curvature will behave like R ~ ce ^ where h and c can be defined from the 
dominant term in l|A.2|) : 

e M « Na N R N ~ l « NaN^'h^- 1 ^. (A.5) 

Here the requirement /' > allows for the following sign combinations of the coefficients a^ and the curvature 
asymptotics R{4> — » oo): 

N = 21 : sign [ojv] = sign [R(<f> — ► oo)] 

TV = 2Z + 1 : aw > 0, sign [£(<£ -> oo)] = ±1. (A.6) 

The other combinations, N = 21 : sign [a if] = —sign [R{(f> — > oo)], A^ = 2/+1 : a^ < 0, sign [R(<j> — > oo)] = ±1, 
would necessarily correspond to the /' < sector, so that the complete consideration should be performed in 
terms of the extended conformal transformation technique of Re f. I15j . Such a consideration is out of the scope 
of the present paper and we restrict our attention to the cases l|A.6|l . The coefficients h and c are then easily 

derived as h = A/(N — 1) and c = sign (o/v) \Ncin\ jy_1 . Plugging this into l|A.3|) one obtains 



(iV-1), , !_ 

. ,.. . .... , ., V J- I AT" "I • 

and that the exponent 



U(<f> -» +oo) » sign ( aAr )Al^_^ ITVajvP^ e~— A *e— A ^ (A.7) 



^- 27V A (A.. 



00 


for D > 2N , 


1 


for D = 2N, 





for D <2N. 



(D-2)(N-1) 
changes its sign at the critical dimension D = 2N: 

U{(t> ^+oo)^ sign (a N )t-^l \Na N \-^ .■ ■{ 1 l«.r /) 2V. (A.!)') 

This critical dimension D = 2N is independent of the concrete coefficie nt aj y and is only defined by the 
degree deg^(/) of the scalar curvature polynomial /. From the asymptotics (|A.9J| we read off that in the high 
curvature limit </> — > +oo, within our oversimplified classical framework, the potential U(<f>) of the considered 
toy-model shows asymptotical freedom for subcritical dimensions D < 2N, a stable behavior for <xn > 0, 
D > 2N and a catastrophic instability for a^r < 0, D > 2N. We note that this general behavior suggests a way 
how to cure a pathological (catastrophic) behavior of polynomial R x — nonlinear theories in a fixed dimension 
D > 2N\. By including higher order corrections up to order N 2 > D/2 the theory gets shifted into the non- 
pathological sector with asymptotical freedom. More generally, one is even led to conjecture that the partially 
pathological behavior of models in supercritical dimensions could be an artifact of a polynomial truncation of 
an (presently unknown) underlying non-polynomial f(R) structure at high curvatures — which probably will 
find its resolution in a strong coupling regime of M— theory or in loop quantum gravity. 



B Graphical visualizations 

Following section |3 we consider the model with one internal space and critical dimension D = 8 (as usual, for 
the external spacetime D = 4). Then, the effective potential l|2.20[l reads: 



U eff (J3\<t>)=e 



— v-40 1 



±R ie -^ + U{<t>) 



(B.l) 



To draw this effective potential, we define R\ via U(<fio) in Eq. I|2.25[l . In its turn, U(<po) is defined in Eq 

H3.24JI . In Figs. 1,2 the generic form of the U e f 
-0.1 from the stability regions lpT2l)j) and (|?07|) 



(I2.:M|) where R(<f>o) = R^l and /' [R{o)i) can be found from 1J3.24JI . In Figs. 1,2 the generic form of the U ef f 



is illustrated by a model with parameters a = 1, 7 = 1, Ag 




Figure 1: Typical contour plot of the effective potential U e ff(f3 1 ^ 
1,7 = 1, A 8 = —0.1. U e ff reaches the global minimum at (/3 1 = 0, ( 



;iven in Eq. IjB.lfl with parameters a 
-2.45). 




Figure 2: Typical form of the effective potential U e ff(f3 1 ,c/)) given in Eq. (|B.1|) with parameters a = 1,7 
l,Ag = -0.1. 
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